Solid liquid phase changes with different densities 
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Abstract. In this paper we present a new thermodynamically consistent phase tran- 
sition model describing the evolution of a liquid substance, e.g., water, in a rigid 
container Q when we freeze the container. Since the density Q2 of ice with volume 
fraction /3 2 , is lower than the density Q\ of water with volume fraction f3\ , experi- 
ments - for instance the freezing of a glass bottle filled with water - show that the 
water pressure increases up to the rupture of the bottle. When the container is not 
impermeable, freezing may produce a non-homogeneous material, for instance water 
ice or sorbet. Here we describe a general class of phase transition processes including 
this example as particular case. Moreover, we study the resulting nonlinear and sin- 
gular PDE system from the analytical viewpoint recovering existence of a global (in 
time) weak solution and also uniqueness for some particular choices of the nonlinear 
functions involved. 

Key words: phase transitions with voids, singular and nonlinear PDE system, global 
existence of solutions. 
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1 Introduction 

In this paper we consider a phase change process occurring in a bounded container 
Q of M 3 . More in particular, we consider a liquid substance (with density Qi) in a 
rigid and impermeable container Q and freeze it. We assume that the solid phase 
has zero velocity, while the liquid phase velocity is the vector Ui , let (3\ be the 
local proportion of the liquid phase, j3 2 the local proportion of the solid phase (with 
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density g 2 < Qi), and § the absolute temperature of the system. In this setting 
(cf. the following Section [2] for further explanations on the model), it is possible to 
observe macroscopically the formation of voids, thus we do not have the condition 
Pi + P2 — 1) but we only know that (3\ + p 2 G [0,1]. In the following Section [2] we 
introduce the model which will lead to the following PDE system we study - coupling 
it with suitable initial and boundary conditions - from the analytical viewpoint in the 
next Sections [3HS1 We search for solutions Pi, @2,p) (p standing for the pressure) 
of the Cauchy problem associated with following PDE system, coupled with suitable 
boundary conditions, in the cylindrical domain Qt := Q x (0,T), being T the final 
time of the process: 

(1.1) 
(1.2) 

tf~ ( *~* c) ] + fy>(A, ft) 9 0. (1.3) 

/ 

Here A represents the Laplacian with respect to the space variables and we have 
assumed that k, £ and $ c are three constants denoting, the compressibility of the 
substance, the latent heat and the critical temperature of the phase transition pro- 
cess, respectively, while R stands for an entropy source. Moreover, v represents 
a non-negative interfacial energy coefficient and tp is a general proper, convex and 
lower-semicontinuous function, taking possibly value +oo in some points. Hence, d(f 
represents here the subdifferential of Convex Analysis (cf., e.g., [20]) and it takes 
possibly into account of the restrictions on the phase field variables Pi and p 2 . In- 
deed, it may be the indicator function of the set of admissible values of Pi and p 2 . 
If K = {(Pi, P2) £ [0,1] I Pi + P 2 G [0,1]}, then a possible choice in our analysis is 
(p(Pi, P2) = Ik(Pi,P2), which is defined to be if (Pi, P 2 ) £ K and +00 otherwise. 
This is the reason why we have an inclusion symbol in (ll.3p . since the dip may possibly 
be a multivalued operator. 

We refer to Section [2] for the detailed description of of our model, but we prefer 
to anticipate few words here on the derivation of equations f ll.lHl.3l) . Equation (II. ip 
follows from the mass balance (cf. also [16]) and the equation of macroscopic motion, 
which is here assumed to be in agreement with the standard Darcy law (cf. [14]). 
Equation (II. 2p represents the rescaled internal energy balance, which is a consequence 
of a new approach developed in [6] (and studied also in [3] [5], [I] [9], and [IE]). Note 
that the presence of the logarithmic contribution in ( 11.211 introduce new difficulties 
in the mathematical analysis of the problem, due to its singular character, but, once 
one has proved that a solution exists in some proper sense, this immediately entails 
the positivity of the temperature, which could otherwise represent a difficult step in 
the proof (cf., e.g., [13]). Inclusion ( 11. 3D stands for the microscopic motion equation 
(cf. pj)] and also [H] and [22] for different approaches) for a proper choice of the 
free-energy functional and of the pseudo-potential of dissipation which, following the 
theory of [20] , accounts for dissipation in the model (cf . also the next Section [2] for 
further details on this topic). 



kd t p + QidtPi + Q 2 d t p 2 - Ap = 0, 
d t (logtf) + ^d t Pi - A(logtf) = R, 
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In all references quoted above (but in [TU]) only phase transition processes in 
which no voids nor overlapping between phases can occur are studied. In that case 
the proofs mainly rely on the fact that (3% + (3% — 1 , hence the mass balance equation 
as well as the effects of the pressure can be neglected and so the system reduces to 
the coupling of equation ( 11. 2ft (with only one phase variable, say \ — Pi ~ fii ) an d 
(11.31) . More in particular, in [B] the model coupling these two equations and accounting 
also for some memory effects in (jl.2p is studied and an existence - of weak solution 

- result is proved in case (p = ip ji , while uniqueness remains still an open problem. 
In [7j and [U [5] uniqueness of solution is established and also some investigation on 
the long-time behaviour of solutions is performed respectively in case v = and 
without any thermal memory for a system in which a different heat flux law in (jl.2p 

- leading to an internal energy balance with d instead of \ogd inside the Laplacian - 
is taken into account. Moreover, in [5] (and in [3] for the case v = and without any 
thermal memory contribution) the long-time behaviour of the solutions-trajectories 
to an extension of the model presented in [6] is studied, while in [8] a more general 
approach to this kind of processes - leading to more general type of nonlinearities in the 
equations - is performed. In particular in [8] the log^-term in ( II. 2p is substituted by 
a more general nonlinear function (maybe also multivalued) including the logarithmic 
function as a particular case. An existence-long-time behaviour of solution result is 
established in [8]. Only in [16] - instead - a first attempt to approach phase-field 
systems including the possibility of having voids is developed. In jT6j indeed we 
study a system coupling the phase equation (11.31) (including also a viscosity term), 
the temperature equation (11.21) . a general mass balance equation, and a quasi-static 
standard stress-strain relation, including viscosity effects but no accelerations. We 
establish in [1Q\ a global well-posedness result for the resulting PDE system coupled 
with suitable initial and boundary conditions. 

Here, under the assumption that the solid phase velocity is and by means of the 
Darcy law (cf. [H] and also (I2.16P below), the mass balance equation can be rewritten 
as ( I3.15p . The main mathematical difficulties stand in the nonlinear and singular 
coupling between ( I3.16P and ( 13.171) . Our approach to this system is to regularize the 
log-term in ( 11.21) by introducing its Lipschitz-continuous approximation. Hence, we 
first solve the regularized problem in Section HJ getting also uniqueness of solutions, 
then, in the last Section we prove our main result (stated in Section [3]) entailing 
global (in time) existence of (weak) solutions to (ll.lHl.3p coupled with suitable initial 
and boundary conditions. The problem of uniqueness for system ( ll.H - ll.3l) is still open 
(cf. also Remark 13.81 below) and this is mainly due to the lack of regularity of the 
■^-component of the solution ($, /3 2 ,p) to (jl.ltil.3p . 

2 The model 

In this Section we detail the derivation of the model leading to the PDE system 
( ll.ltil.3p already introduced in Section CD Hence, let us consider a liquid substance, 
e.g., water, in a rigid and impermeable container Q and freeze the container. Since 
the density g 2 of ice with volume fraction (3 2 is lower than the density £1 of water 
with volume fraction (3\ , it seems impossible that water transforms into ice (or, more 
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in general, into a mixture of ice and water): the water remains liquid even at low 
temperature. Experiments - for instance the freezing of a glass bottle filled with water 
- show that the water pressure increases up to the rupture of the bottle. When the 
container is not impermeable, freezing may produce an non homogeneous material, for 
instance water ice or sorbet. In this setting, it is possible to observe macroscopically 
the presence of voids and thus we do not have the condition ft + ft — 1 , but we can 
only deduce that < ft + ft < 1 . 

We assume the ice has zero velocity, while the water velocity is represented by 
the vector Ui . For the sake of simplicity, we assume the small perturbation assumption 
(cf. p2]). In this setting the mass balance is 

fadtQx + fcdt&i + QidtPi + £>2<9 t ft + Q\{h divUi = 0, 
- QtfcUi ■ n = w, (2.1) 

where vo is the liquid water intake and by n we denote the outward unitary normal 
vector to dVt. We have assumed the substance to be compressible, even in some cases 
this compressibility is very low. Moreover, the equation of motion for the macroscopic 
motion of the water can be written as 

Qid t (ftUi) = divcri + f x , 

am = gl , (2.2) 

where <j\ is the water stress, fi is the action of the exterior on the liquid water, 
for instance the friction of the water on the solid ice phase, and gi is the boundary 
exterior action on the liquid water. 

For the microscopic forces responsible for the ice water phase change, the equa- 
tion of motion is (cf. [To] ) 

Bl Vdivf 5 1 1=0. 



B 2 J \ H 2 

This vector equation is the equation of microscopic motion and it deals with the works 
produced by the microscopic motions occurring in the phase changes. The terms B\ , 
B 2 , H± , and H 2 are new internal forces responsible for phase transitions phenomena. 
In particular Bi represents the energy density per unit of ft , B 2 represents the 
energy density per unit of [3 2 , while if i and H 2 represent densities of energy flux. If 
we assume that no exterior source of work is present, we also have 



Hi 
H 2 



■ n = 0. (2.3) 



The internal energy balance reformulated in terms of the entropy s is (cf. [6] for its 
derivation) 

d t s + div Q = R, (2.4) 



where we have denoted by Q the entropy flux, by R the entropy source, r = R§ being 
the heat source, and by f3 the column vector of components ft, ft. No high order 
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dissipative contribution is involved on the right hand side due to the small perturbation 
assumption. 

The internal forces are split into dissipative forces indexed by d and non dis- 
sipative interior forces indexed by nd . The non dissipative forces are defined through 
the water free energy m which takes into account the internal constraint on the state 
quantities, for instance 

Wi, fo)eK = {( 7l , 72 ) G [0, 1] | 7i + 72 e [0, 1]}. 

We choose the following form for the free energy: 



- co^iogtf - 1) - ^(tf - C ) + MA, AO 

+ \ (|V/3i| 2 + |V/? 2 | 2 ) + f3M gi ) + #2/2(02) 



(2.5) 



where $ c is the critical temperature of the system, I denotes the latent heat of the 
melting solidification process, c the positive specific heat, v a positive interfacial 
energy coefficient, and Ik stands for the indicator function of the set K , it takes 
value on K and +oo otherwise. Because we have assumed the substance to be 
compressible, the free energy depends on the densities. Indeed, fx and / 2 account, 
respectively, for the dependence on the densities of the free energy related to the liquid 
and the solid phases. Thus we have 



dm 



s = —- 



c log 1? + 



and we define the non dissipative internal forces as: 



(2.6) 





/ 


\ 


/ B nd \ 


dBx 


V Bf J ~ 


dm 








J 



nd 



HI 



PI 



nd 



pnd 
r 2 



( 



dm 



d(V(3x) 

dm 



\ 



v 



\ d{V(3 2 ) ) 

I dm \ 

dgx 

dm 
\~dV2' 



V(3x 
V/5 2 



dpi 
3 dh 



+ dI K (Px,P2) 



(2.7) 



(2.8) 



(2.9) 



with an abuse of notation for the derivatives dm /d(3 and where dlx represents the 
sub differential of Ik in the sense of Convex Analysis (cf., e.g., [TU] or [2D])- Because 
the free energy does not depend on deformations, we have 



a 



nd 



0. 



(2.10) 
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The dissipative forces are defined with pseudopotential of dissipation $ which takes 
into account the internal constraints related to the velocities and ensures that the sec- 
ond law of thermodynamics is satisfied. The mass balance represents such a constraint 
(cf. also [16J for further comments on this choice). For instance, we may choose 



= I (PxdtQx + (3 2 d t g 2 + QidtPi + Q2dA + Q1P1 divXJx) 



+ ^(Wi\ 2 + \d t p 2 \ 2 ) + ^m 2 , 

where Iq denotes the indicator function of the origin of R, // the phase change vis- 
cosity and A > denotes the heat conductivity of the system, while the operator 
D represents the linearized symmetric strain tensor £)y(U) := (Uij + Uj j i)/2, where 
i,j = 1,2,3 and with the commas we mean space derivatives. The dissipative forces 
can be written as: 



A 



{Pf, P d , Bf, B d , Hi Hi at Q d ) 

e <9$ (dtQt, d t g 2 , d t Pi, d t [3 2 , V(<9 t /?i) 

where <9$ represents the sub differential of $ and it results 

3$ 



V(d t /3 2 ), D(Ux), W) 



Q = Q C 



( 



( pf\ 

JDd 
r 2 

Bf 
B d 2 

V °i I 



0, H% 



<9$ \ 



-AV(logtf), 



d{d t Qx) 
d(d t g 2 ) 



d(Vd t p 2 ) 










0. 



(2.11) 
(2.12) 



dl (0) 



0i 
ft 
Qi 



with an abuse of notation for the derivatives d§/d(d t X) and where 1 stands for the 
unit tensor. By denoting — p an element of dlo , the last relations can be written as: 



2; 



Pf = -p/3 u P d 
B d = fidtfc - pqi, B d = nd t (3 2 - pg 2 , 

of = -vexfcl- 



(2.13) 
(2.14) 
(2.15) 



The quantity pQ\(3i turns out to be the pressure in the liquid phase and using consti- 
tutive law (12.15ft in the equation for the macroscopic motion (12.21) we get 
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We may assume that the acceleration is negligible and that the volume exterior force 
fi results mainly from friction on the solid phase and that this force is proportional 
to the relative velocity through the relation 

m 

This is mainly due to the fact that the force has to be when there is no water. 

Remark 2.1. The density of force fi with respect to the actual volume of water flidQ 
is 



fi e?fl 



Pi m 

It is proportional to f3\ and to the liquid water momentum £>i/3iUi as one may expect. 
With this choice, from (12. 2\\ we get 

m 

which (within the small perturbation theory (cf. [17J)) can be rewritten as: 

Ql (3 1 \J 1 = -mVp, (2.16) 

where the parameter m denotes the mobility of the water in agreement with the Darcy 
law (cf. [i4j). 

Since the densities Q\ and q 2 are not related to a particular internal force, their 
constitutive laws are (cf. [El (3.25), p. 11]) 

pud + pd = 0> pud + pd = Q) 

which, with (12.91) and (12.131) . can be written as: 

a P = >x P = - 217 

OQi OQ 2 

Within the small perturbation assumption, we have 

p = ki(Qi - Q r i),P= k 2 (62 - £2); ( 2 - 18 ) 

where the ki (i = 1, 2) are the compressibilities of the liquid and solid phases, respec- 
tively, and the g[ (i = 1, 2) are reference densities. Let us note that the densities, 
which depend on the pressure p through (12.181) . may also depend on the temperature 
via a thermal expansion coefficient, giving p = ki(gi — q[) + kf 1 ^ — d c ) in place of 
(12.18!) . However, for the time being, we do not introduce this dependence. 

The predictive theory equations within the small perturbation assumption result 
from the mass balance (12.11) with (12.181) . the equation of motion (I2.2H2.3I) . the entropy 
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balance (12. 4p and the constitutive laws (I2.6H2.14]) with relation (12.161) . Hence, they 
can be written as follows 

kd tP + Ql d t 1 + Q 2 d t 2 - mAp = 0, (2.19) 

coa t (logtf) + jd t (5 x - AA(logtf) = R, (2.20) 

M ( ** ) -A )~P ( V ( -* ( * " * J 1 + 3 0, (2.2!) 



d t fo J \A/3 2 J F \g 2 J ' \ 

where k = I— — h — J is the compressibility of the mixture, being /3° + the 

reference value of the material volume fraction. This set of equation is completed by 
the following boundary conditions on dQ 

= 0, -p- = 0, 2.22 
an an 

A— + oi${-& — ^ext) — 0, (2.23) 
an 

where a$ is the thermal conductivity of the container boundary. Concerning the 
hydraulic boundary conditions, we may choose 

dv 

m^- + a p (p - Pext ) = 0, (2.24) 

where a p = if the boundary is watertight, a p = oo if the boundary is connected to a 
water supply with pressure p ext (in this case the boundary condition is p = p ex t)- This 
condition means that the flow is proportional to the difference of pressure p — p ext , in 
case the container is permeable. In case, instead, of a metallic recipient, the condition 
accounts for the deformation of the metal due to the pressure. 

Remark 2.2. Due to the small perturbation assumption, densities q% and g 2 in 
equations (12.19)) and 112. 2 lp are constant. 

Remark 2.3. More sophisticated boundary condition result from more sophisticated 
physical boundary properties. The boundary of the container can be semi-permeable, 
in this case we have 

m- OL v {p-p ext Y = 0, 

an 

where (p — p e xt)~ = sup(0, — (p — p ex t)) is the negative part of p — p ex t ■ This boundary 
condition means that when the pressure is lower than the exterior pressure p ext , water 
flows inside the container but when the pressure is larger than the exterior pressure, 
no water flows outside. 

The pressure may also be controlled on the boundary by the following relation 

dp 

m— + dl-(p - Pext) 3 0. 
an 

Water flows outside the container in order to maintain the pressure lower than the 
outside pressure p ext which may be the atmospheric pressure. 
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The initial conditions in f2 are 

P x (x,0) = ft(x), &(x,0) = /? 2 °(x), 0(x,O)=& o (x), p(x,0)= Po (x). 

This set of partial differential equations allow to compute pressure p(x, t) , liquid water 
and ice volume fractions fli(x,t), /3 2 (x,t) and temperature $(x,t) depending on the 
external actions resulting from the exterior pressure p ex t{z, t) and exterior temperature 
$ ext (x, t) and rate of heat production $ c R(x, t) very often equal to . In an engineering 
situation indeed the governing action is the exterior temperature cooling and heating 
the container. 



2.1 Some examples 

1. Freezing water in an impermeable container. Suppose to have = 1, fi\ = 
0, $0 > and that the boundary condition for the pressure is ( 12.24ft with a p = 0, i.e.: 

£ = 0. 
an 

Hence, the container is completely cooled. In order to look for closed form solutions we 
assume all the quantities are homogeneous. Assuming the temperature to be known, 
we focus on equation (12.211) which describes the phase change. Because of homogeneity, 
it is 

"(2£)-'(S) + ("* < o"'" ) ) + ^« att 

When the Celsius temperature $ c is positive, the solution of this equation is fli = 1 
and /3 2 = 0. When the Celsius temperature becomes negative, the phase change which 
is expected does not occur. Indeed due to the mass balance, to values Pi = 1 and 
@2 = and to condition /3i + /3 2 < 1 , we should have, neglecting the effect of the 
compressibility, the following relations: 

Q2d t (3 2 + QidtPx = 0, dtfa < 0, dtfo + <9 t /3i < 0. 
Note that the two first conditions imply 

dtfc + <9 t A = - dtfa > 0, 

which contradicts the third condition. Let us check that the water remains liquid 
(dtfli = 0) even if Celsius temperature is negative. We must have 

with 

dI K (l,0)=(^Y P>0,Q<P, 
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or 

PQi + ~ A) > 0, 
0<p(f?i-f? 2 ) + ^(^-^ c ). 

Since Qi — g 2 > 0, these conditions may be satisfied. In this case, the minimum value 
of the pressure is 

P=— f ^y(tf-tfc). 

[Ql - Q2)V C 

Thus, when cooling the container, water remains liquid and pressure increases. This 
is in agreement with experiment: a glass container filled of water explodes due to the 
pressure increase when freezing. 

2. Freezing of a water emulsion in an impermeable container. Suppose to 
have < (3® < 1, (3 2 = 0, $0 > 0, hence, the emulsion (a mixture of voids and water) 
is cooled. Assuming an homogeneous evolution, it may be seen that the water freezes 
till either /3i = 0, (in this case (3® < g 2 j g\ and the cooling results in a mixture of 
ice and voids), or @i + (3 2 — 1> ( m this case /3° > g 2 j Qi and the cooling results in a 
mixture of ice and liquid water, (3\ = (gift® — Q2)/{.Q\~ 82), P2 = (Qi~ Qifii) I (Qi — 82))- 
This mixture is a water ice or sorbet. 

3. The ice in a glacier. Assume to have an homogeneous mixture of ice and water 
without voids at an equilibrium, 

ft +& = 1, 8^ = 8^2 = 0, 

then, equations fl2.19ti2.2T]) give 

-(£W"* ( o-' J ) + (? 

with P > . Hence, we get 

-pG?i-f? 2 )-^-(tf-tf c ) = 0, P >o. 

In this situation the phase change temperature depends on the pressure. It is known 
that this is the case for the ice-water phase change. Because g± > g 2 the phase change 
temperature decreases when the pressure increases. A consequence of this property is 
that the ice melts at the bottom of a glacier and lubricates the rock ice contact surface 
allowing the downhill motion of the glacier. 



3 Variational formulation and main results 



In the following part of the paper we want to deal with a suitable generalization of the 
PDE system ( l2.19ti2~!Tn) introduced in the previous Section [21 Hence, we first detail in 
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this section the assumptions on the data and on the general proper, convex, and lower 
semicontinuous mapping ip : R 2 — > [0, +00] , which will generalize the role played by 
the indicator function Ik in the inclusion (12.21j) introduced in the previous Section El 
and then we state a suitable variational formulation of the generalized PDE system 
obtained and our main results. 

In order to do that, we first introduce the Hilbert triplet (V, H, V) where 

H := L 2 (VL) and V := H 1 ^), (3.1) 

and Q C R 3 is a bounded and connected domain with Lipschitz continuous boundary 
T := dfl. Let T be a positive final time of the process and denote by Qt ■— Q x (0, t) 
and S; := T x (0, t), t G (0, T]. Then, we identify, as usual, H (which stands either 
for the space L 2 {Q) or for (L 2 (f2)) 3 or for (L 2 (fi);R 2 )) with its dual space H' , so 
that V ^ H "—>■ V with dense and continuous embeddings. Moreover, we denote by 
|| ■ \\x the norm in some space X , by (*, •) the scalar product in H , and by (-, •) the 
duality pairing between V and V . Hence, for any £ G V , set 

C n :=p<C,l>, (3-2) 
V := {C E V : (n = 0}, V:=VnV. (3.3) 

The above notation V is just suggested for the sake of convenience; indeed, we mainly 
see V, V as (closed) subspaces of V, V , inheriting their norms, rather than as a 
couple of spaces in duality. 

^From now on, for simplicity and without any loss of generality, we suppose 
that the coefficients in (I2.19l - [2~?i1) k = m = co = \ = fi = £/i9 c = l. Next, in order to 
give a precise formulation of our problem f l2.191 - [2~2TT) . we define here the realization 
of the Laplace operator with Neumann homogeneous boundary conditions, that is the 
operator 

B : V -»• V , (Bu, v)= [ Vu-Vv u,veV. (3.4) 

Jn 

Clearly, B maps V onto V and its restriction to V is an isomorphism of V onto V . 

Define W := {w G H 2 (Q) : d n w = on T}, where d n is the derivative with 
respect to the outward normal derivative to T and make the following assumptions on 
the data. 

Hypothesis 3.1. We assume that Q\ > Q2 > 0, v > 0, and suppose that 

: M. 2 — > [0, +00] is proper, convex, lower semicontinuous and (p(0) = 0, (3.5) 



V(<p) bounded iiv = 0, (3.6) 

7(7") := exp(r) for r G R, (3-7) 

iJoGl 1 ^), tf >0 a.e. infi, w := 7 _1 (^o) € H, (3.8) 

/3 = (/O2) e vfleV, Po eV, (3.9) 

R G L 2 (Q T ) fl L 1 (0, T; L°°(Q)), n G L°°(S T ). (3.10) 
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Then, we introduce the functions 1Z G L 2 (0,T; V) such that 

(K(t),v) = I R(t)v+ I U(t)v ldn veV, for a.e. t G [0, T\. (3.11) 
Jn Jan 

Moreover, we denote by dip the sub differential in M 2 x JR 2 of the convex analysis 
and note that dip is maximal monotone and that d(p(0) 3 (see, e.g., [2] and [TU] 
for the general theory). The same symbol dip will be used for the maximal monotone 
operators induced on L 2 -spaces. 

Remark 3.2. Here in our analysis we have chosen to treat the Neumann homogeneous 
boundary conditions on p, that is, we consider (12.241) with a p = 0. It is not difficult, 
however to treat the boundary conditions f)2.24p with a p > with the same tech- 
niques used here. Regarding, instead, the absolute temperature i? , we work here with 
Neumann non-homogeneous boundary conditions on log i? , which means we consider 
(12.23!) m case «tf > and f} ext = . Other type of boundary conditions like Dirich- 
let non-homogeneous boundary conditions on log d could also be taken into account 
without any difficulties in our analysis. Regarding instead the more sophisticated 
boundary conditions on p detailed in Remark 12.31 a more careful analysis should be 
performed in order to include them in our results, but we do not want to face this 
problem in this paper. 

Then, we are ready to introduce the variational formulation of our problem as 
follows. 

Problem (P). Find (w, f3i,f3 2 ) and ($, p) with the regularities 

w G -f^ 1 (0, T; V') fl L 2 (0, T; V), $ G L 5/3 (Q T ) (3.12) 

(3 h /3 2 G H\0,T;H), u/3 u u(3 2 G L°°(0, T; V) fl L 5//3 (0, T; W 2,b ^{Vt)), (3.13) 
P=(j3 1 ,p 2 )eT>(<p) a.e. in Q T) 

£i, 6 e LV\Q T ), p G C°([0, T}; V) n L 2 (0, T- W) n ^(0, T- H), (3.14) 
satisfying 

d t p + Qidtfii + Q 2 d$i + Bp = a.e. in Q T , (3.15) 

d t w + dtfc + Bw = TZ inV a.e. in [0,T], (3.16) 

^-P +? - p fe)K""o" C ) a ' em *' (317) 

= 7 ( w ), I = (a, 6) G 9^09) a.e. in Q T , (3.18) 

and such that 



w(0) = wo, p(0) = po a.e. in fl, 
/3(0) = (/5 1 (0),/5 2 (0)) = /3 a.e. in a 



(3.19) 
(3.20) 
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Remark 3.3. Notice that testing formally equation (13.151) by 1 one can immediately 
deduce that the "total mass" of p + f3 e , where f3 e := (qiPi, 02P2) is conserved during 
the time interval [0, T] . Indeed, we have 

(p + QxPi + Q2P2) (t) = f(p+ Q1P1 + Q2P2) (0) (3.21) 

Jn 

for all t G [0, T] . The physical meaning of this relation is that the mass of water 
either liquid or solid contained in the domain Q is constant, due to the watertight 
boundary condition we have chosen (cf. ( 12.24ft with a p = ). Finally, let us note that, 
by interpolation, from (I3.12p it follows that w G C°([0, T];H) and so all the three 
boundary conditions in h3.1^\3~W^i hold in H . 

Remark 3.4. Obviously, an admissible convex function satisfying assumption ( 13. 5ft 
is, e.g., if = Ik , being K the convex set of IR 2 defined in the Introduction. 

Remark 3.5. Let us underline here (cf. also |6[ Remark 3.1]) some regularities prop- 
erties of possible solutions to Problem (P). First of all, let us recall the Gagliardo- 
Niremberg inequality in 3D, that is 

\\v\\ L P { n) <c|M|£|M|^ Q VveV, (3.22) 

for ^ = § + ^ip and for some positive constant c . Hence, if we are able to prove that 
e l°°(0, T; H) D L 2 (0, T; V) , then it follows that 

§ G L 5/3 (Q). 

Remark 3.6. Note that in our analysis we could also consider a more general form 
of our equation ( 13.1 7ft ; we can treat the case in which an antimonotone, but smooth, 
contribution is added to the monotone one given by dip . Indeed, we can cover with 
our analysis the case in which (13.171) is substituted by the more general inclusion: 

dtP + v + d V {(3) + & '03) - p Q 9 ft ~ ^ a.e. in Q T , (3.23) 

where a G C 1 ' 1 (M 2 ) is a, possibly antimonotone, contribution coming from the possibly 
non-convex part c(/3) in the free energy functional (12. 5p . The example we have in 
mind is the one given by the "classical" double-well potential. 

The aim of the following Section 15.11 is to prove one of our main results, that is 
the following global existence theorem. 

Theorem 3.7. Let Hypothesis 13. II hold true and let T be a positive final time. Then 
Problem (P) has at least a solution on the whole time interval [0, T] . 

Remark 3.8. The problem of finding also a uniqueness result for Problem (P) is 
still open and is strictly related to the possibility of finding more regularity on the 
$ -component of the solution to Problem (P). The reader can refer also to IE$ and 
0/ for further comments on this topic in case the voids are not admissible and only 
the two equations (13.161) and (I3.17P (with p = and ip = Iq, being C = {(Pi, P2) G 
[0, 1] I fli + @2 = 1}) are coupled. 
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The proof of Theorem 13.71 is based on the following scheme, which follows the 
idea of [BJ: we first solve an approximating problem in which 7 is substituted by a 
Lipschitz-continuous function, then we make a priori estimates (independent of the 
approximation parameter), which will allow us to pass to the limit by means of com- 
pactness and monotonicity arguments. In view of these considerations, let us state 
here a preliminary result (whose proof will be given in Section H]). 

Consider the following assumptions on the data 

Hypothesis 3.9. Let A3.5U3.6}) and A3.y\) hold true and assume moreover that: 

5 : = 7" 1 : R -f R 2 is maximal monotone operator 

with Lipschitz continuous inverse graph 7, (3.24) 

w eH, Tie L 2 (0,T;V), (3.25) 

where 1Z is defined in A3. 

Theorem 3.10. Suppose now that Hypothesis 13.91 holds true. Then Problem (P) 
has at least a solution (w, Pi, P2) , £,i,£,2, V) satisfying the following regularity 
properties: 

w,$e C°([0,T};H)nL 2 (0,T;V), (3.26) 
vPi, vp 2 e L 2 (0,T;W), ^1, 6 e L 2 (0,T;H). (3.27) 

Moreover, the components w , Pi, p 2 , p, and § of such a solution are uniquely deter- 
mined. 

Remark 3.11. Observe that the main advantage of taking the entropy balance equa- 
tion (I3.16P instead of the internal energy balance equation is that once one has solved 
the problem in some sense and has found the temperature d := 7(w) , it is automat- 
ically positive because it stands in the image of the function 7 (cf. ( 13.71) ). Indeed 
in many cases it is difficult to deduce this fact only from the internal energy balance 
equation (cf, e.g., \13] in order to see one example of these difficulties) . Let us note 
that within the small perturbations assumption the entropy balance and the classical 
heat equation are equivalent in mechanical terms (cf. 

Moreover, in the same framework as above, let us recall the following inequality, 
holding in Q introduced above (cf. (13. 2p ): 

\\v-v n \\ H < c\\Vv\\ H \/v G V, (3.28) 

for some constants c > 0, depending only on Q. Inequality (I3.28P is one form of 
the standard Poincare-Wirtinger inequality. Finally, we will make use of the following 
elementary inequality: 



ab < na 2 + —b 2 , Va,beR V77 > 0. 
4?7 



(3.29) 
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4 Proof of Theorem 13.10 

In this section we assume Hypothesis 13.91 and we prove Theorem 13.101 Here we are 
going to prove existence and uniqueness of solutions to Problem (P) £ in the time 
interval [0, T] , using a fixed point technique of contractive type. 

In order to develop the proof in a rigorous way, we introduce here the Lipschitz 
continuous Yosida-Moreau approximation dip £ = (ip £ )' (cf. [TQJ, Prop. 2.11, p. 39]) of 
dip, we first solve the Problem (P) with d(p £ instead of dip (call it (P) £ in this 
section) and then we perform a priori estimates independent of e and pass to the 
limit as e \ recovering a solution to Problem (P). Moreover, we denote by c the 
positive constants (maybe different from line to line), depending on the data of the 
problem, but not on e . Then we are ready to introduce the variational formulation of 



our approximating Problem (P) £ as follows. 

Problem (P) £ . Find (u> £ , /3f,/3f) and ($ £ , £f,£f, p £ ) with the regularities 

w £ } tf e G H X %T\V') nL 2 (0,T;V), (4.1) 
(3 e = ((3l(3 £ 2 ) G (H\0,T;H)) 2 , v(5\, vft G L°°(0, T; V) fl L 2 (0, T; W), (4.2) 

e L 2 (Q T ), p £ e C°([0,T] ] V)nL 2 (0,T;W)nH\0,T ] H), (4.3) 

satisfying 

d tP £ + Q X d t p{ + g 2 d t (3 £ + Bp £ = a.e. in Q T , (4.4) 

d t w £ + d t (3l + Bw £ = TZ in V, a.e. in [0, T], (4.5) 

^ + " (el) + « E " {%) = (' o " c ) a - in < 46 ) 

r = 7 K), r=(S,«5) = (^)W a.e. inQ T , (4.7) 
and such that 

w £ (0) = w , p e {0) = po a.e. in Q, (4.8) 

f3 £ (0) = (/3f(0),/3|(0)) = /3 a.e. in fi. (4.9) 



Now we start here proving existence of solutions to (P) £ , by means of a con- 
traction argument. 

Existence of solutions to (P) e , e G (0, f ] . Let us take t G [0,T] (we will choose 
it later) and call 

X:={(r,8)e(H 1 (0,P,H)) 2 : (r,s) eV((p)}. 

First fix (fli, fa) E X in the equations (14.41) and f)4.5p . then, by well-known re- 
sults (cf. [H Lemma 6.3] and also (13.101) ). we find also a unique w := 7i(/?i,/? 2 ) G 
H\Q,T;H) n C°([0,T];1/) n L 2 (0,T; H 2 (tt)) solution of g3D and p := T 2 (p u 2 ) G 
J ff 1 (0,T; J ?J)nC ([0,r];1/)nL 2 (0,T;H/) solution to (03} • Then, if we take this values 
of w and p in (j4.6p in place of w e , p e , again by standard results, we can find the 
solution (3 = (f3x,f3 2 ) to (Oil . 
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In this way, we have defined an operator T : X — > X such that (/^i,/^) =: 
T(Pi,P2)- What we have to do now is to prove that T is a contraction mapping on 
X for a sufficiently small £ G [0, T] . In order to prove that T is contractive, let us 
proceed by steps and forget of the apices e. 

First step. Let e Af (z = 1,2), = 7^(/3/,/^), = T^,^), and 

= TiPiJi). Then, writing two times @3) with (Pi,^) (« = 1,2), making 
the difference, testing the resulting equation with (w 1 — w 2 ) , and integrating on (0, t) 
with t G [0, i], we get the following inequality 

Ik 1 ^) - ^ 2 Wlli + Ik 1 - ^liw) 

< Ci^lK^i 1 )* - + IK^ 1 )* - (/3 r 2 2 )*||| 2(0 ,t;^)), (4.io) 
for some positive constant C\ independent of t. 

Second step. We write now (I4.4p with (/?Y, /3 2 l ), make the difference, test the 
resulting equation with p 1 — p 2 , and integrate on (0, t) with t G [0, t\ , then we get the 
following inequality 

IKP 1 - P 2 )(t)\\ 2 H +\\P 1 -P 2 \\h { 0,t;v ) 

< C 2 (\m\ - {Pi)t\\l^ t]H) + IK^ 1 )* - (/3 2 2 )*||! 2(0it;H) ), (4.11) 
for some positive constant C 2 independent of t. 

Third step. Let us take $ l = ^(w 1 ) (i = 1, 2) and write equation (14.61) for -d 1 and u l , 
make the difference between the two equations written for i — 1 and % = 2, and test 
the resulting vectorial equation by the vector ((/?i)t — (/3 2 )*, (Z^)* ~~ (/3f)t)- Summing 
up the two lines and integrating on (0, t) with t G [0, t] , we have 



E IK' 9 ))* - (z^ii W) + *E n v (# - #)(*)ii 

i=i i=i 



< / |(p 1 -p 2 )(ft((/3 1 1 ) 4 -(/3 1 2 W + fe((/3 2 1 )*-(/3 2 2 )0)| 

< 5 E IK/ })* - Who**) + - ^ 2 n W) + 21b 1 - p 2 ii W)- 

Moreover, using the assumption (I3.24p on 7, we get the following inequality 

2 2 

2 



i E HWft - Whoa*) + ^EllV(/3) - #)(t)|& 
J i=i i=i 

< 2* (llw 1 -wYco m m + lb 1 -P 2 |lco ([ o lt ];if)) • (4-12) 
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Fourth step. Summing up the two inequalities (14.111) and f)4.12p . and using (14.101) . 

we get 

2 2 

E " OS?)*Hi>(0,«;fl) + H V (# " /?)(*)Hh 

i=i i=i 

2 

i=i 

for some positive constant C3 independent of t. Hence, choosing t sufficiently small 
(this is our t), we recover the contractive property of T . Moreover, applying the 
Banach fixed point theorem to T, we get a unique solution for the Problem (P) £ on 
the time interval [0, i\ . Then, it is not difficult to check that 

~~ 1 — 1 & ~~ 2 — ^ k (^-^3-^"')^' ~~ 1 — 1 ~ 2 ~ 2 

,/3 2 ) -T{&\ ,@2 ) \(W(0,T;H))2 < ,fc ) - (A ,/?2 ) | (ffl (0,T;H)) 2 

for any power T fc of T. Hence, one can find k such that {C^T) k < k\ and consequently 
the corresponding T k is a contraction and admits a unique fixed point, which is the 
searched solution to Problem (P) e . on [0, T] . 

Estimates and passage to the limit. We perform now uniform in e estimates on 
the solution to Problem (P) e , which allow us to pass to the limit as e \ 0. 

First estimate. Test (14 .4ft by p, (14.51) by d = 7(10), and (14.61) by d t {3. Then, sum 
up the result of the first two tests with the sum of the two components of the vectorial 
equation found by the third test. This, after an integration in time on (0,t), with 
t G [0, T] , thanks to a cancellation of two integrals, leads to 



2 

2 

v 
2 



db(*)lllr+l|Vp||! W) + / l(w(t)) + c, / / |W| 2 + E Will 

JJ! JO Jn - j 

£l|V#(t)|& + / ^W))<|lbo||^+ / TW+ /V*> 
• =1 Jn z Jn Jo 

+ ^£||V&(0)|& + / ^(/3 )+ f I (4.14) 
z - =1 Jn Jo Jn 

Now, using (13.91) . (I3.241 - [3~25|) . and Schwarz inequality with (I3.29p . we get: 



\\p(t)\\ 2 H + \\mh { 0,t;H) + / Ht)+ / / |W( S )| 2 d S + EWill L 2 (0.t;H) 

Jn Jo Jn . j 

+ II V^(t)||^ + / tp%(3(t)) < c + cJi?||! W0 + r}\m\h( ,t;V), 
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holding true for all 77 > and for some positive c v . Finally, using (13.281) and (13.251) . 
and choosing 77 sufficiently small, we get: 

2 

\\p(t)\\H + l|Vp||i W ) + \\nh(0,t;V) + E IWiH 

i=i 

2 

+ ^||V/^)||^ + / ^W))<c (4.15) 
i=i ^ 

Then, testing (14.51) with w and using (I4.15p . it is a standard matter to deduce the 
following bound 

IMU^(0/r;H)nL 2 (0,T;V)nH 1 (0,T;V'') < c ; (4.16) 

from which, using assumptions (I3.24[) and ( 13.91) . we also deduce 

II^IU^o/r^nL^o^vOnH^ci.T^') + lbl|H 1 (o,T;H)nL° c (o,T ; y)nL 2 (o,T;vy) < c - (4-17) 

Second estimate. Test now (14. 6 p by £ e . Note that the term 

(V(3%t), vr(*)) = (V/3 £ (t), V(^)'(/3 £ (t))) 

is non-negative due to the monotonicity of (y? £ )'. Hence, due to estimate (14.151) and 
(I4.17p . we obtain 

U £ i\\Ll(0,T;H) <C, 2 = 1,2, (4.18) 

and also 

\WPi\\L2{0,T;W) < c, z = 1, 2, (4.19) 

by comparison in (13. 17ft and by standard regularity results for elliptic equations. 

Passage to the limit. As we have just mentioned, we want to conclude the proof 
of Theorem 13.101 passing to the limit in Problem (P) £ as e \ using the previous 
uniform (in e ) estimates on its solution and exploiting some compactness-monotonicity 
arguments. Let us list before the weak or weak-star convergence coming directly 
from the previous estimates and well-known weak-compactness results. Note that the 
following convergences hold only up to a subsequence of e tending to (let us say 
£ k \ 0). We denote it again with e only for simplicity of notation. From the estimates 
(SinHEZD, (HTT8H4TT91) . we deduce that there exist (w,d,p,Q, j = 1,2, such that 



w £ — > w 


weakly star in H\0, T; V) n L 2 (0, T; 1/), 


(4.20) 




weakly star in H l (0, T; V) H L 2 (0, T; 1/), 


(4.21) 


% - ft 


weakly in ^(0, T; #), {j = 1, 2), 


(4.22) 




weakly star in L°°(0, T; 1/) n L 2 (0, T; W), (j = 1, 2), 


(4.23) 


p £ — > p 


weakly star in H\0, T; H) n L°°(0, T; V) fl L 2 (0, T; W), 


(4.24) 




weakly in L 2 (Q T ), = 1,2). 


(4.25) 
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Moreover, employing the Aubin- Lions lemma (cf. [THl P- 58]) and Cor. 8, p. 89], 
we also get 

w £ — > w strongly in f 2 (0,T; if), and hence a.e. in Q T , (4.26) 
■& £ — > # strongly in f 2 (0,T; if), and hence a.e. in Qt, (4.27) 

/3J->/3tj strongly in C°([0,T]; ff) H f 2 (0, T;V) if z/ > (j = l,2), (4.28) 
p £ -> p strongly in C°([0, T] ; if ) fl f 2 (0, T; 1/) and hence a.e. in Q T . 

Note that convergences (14.2 ip . (14.261) . and relation (15.101) imply immediately (cf. [TBI 
Lemme 1.3, p. 12]) the convergence 

$ e — > d = j(w) a.e. in Qt- 

Finally, it remains to prove the identification of the maximal monotone graph dip, i.e. 

£ G d(p(0) a.e. in Q T , 

with £ = ^2) and ^ (j = 1, 2) that are the weak limits defined in ( 14.25}) . In order 
to do that we should verify that 

limsup / T (r,/3 £ )< / T (l,/3). (4.29) 

e\0 Jo JO 

The proof here is splitted into two parts. 



Case v > 0. From ( 14.28}) and ( 14.251) we immediately deduce that ( 15.321) is verified. 

Case v = 0. In this case we can prove that (3 e is a Cauchy sequence in C°([0, T}; H) . 
Indeed one can take the differences of equations (13.171) written for two different indices 
e and e' and test it by the difference vector (f3 e — (3 £ ) . Take the differences of the 
two equations (13.151) integrated in time and test it by p e — p £ . Summing up the two 
resulting equations and integrating over (0, t), t G (0, T) , we get: 

\\p e -p e )\U , t -,H) + h(f3 £ -P £, )(t)\\ 2 H< [ iP-tfm-Pti 

£ JQt 

< fw -^'mum- (3t')(s)\\ H ds. 

Jo 

Applying now the Gronwall lemma ^01 Lemme A. 3], we get 

|| (/3 e - (3 s ')(t)\\ 2 H < c\\F - # £ '\\ LH o,T; H ) - as e, s' \ 0. 

Hence, we have that (3 s — » (3 strongly in C°([0, T]; if), which is sufficient in order to 
prove that (I5.32j) is satisfied. 

Hence, in both cases, all these convergences with the identifications made above 
make us able to pass to the limit (as e \ or at least for a subsequence of it) in 
Problem (P) £ finding a solution to Problem (P). Note that once we prove unique- 
ness, the convergences above will turn out to hold for all the sequence e \ and not 
only up to a subsequence. 
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Uniqueness. Assume that (w 1 , /3 1 ), (i? 1 , p 1 ) and (u> 2 , /3 2 ), (-# 2 , £ 2 , p 2 ) are two 
solutions of Problem (P) corresponding to the same data, and let us call by w : = 
w 1 — w 2 , p := p 1 —p 2 , /3j := /3j — (3 2 , (j = 1,2). Then make the differences from (13 . 151) 
written down for the first solution and H3. 151) integrated in time and written down for 
the second solution, then test the resulting equation by p, test the difference between 
the equations (13.1 7p written in terms of /3 1 and (3 2 by (P 1 ,P 2 )- Finally, integrate in 
time (I3.16P (let us call it 1 * ( I3.16P ) and make the differences from 1 * (13.161) written 
down for the first solution and 1 * (13.161) written down for the second solution, then 
test the resulting equation by w. Then sum up the three resulting equations and 
integrate it over (0,t), use the maximal monotonicity of dip, getting 



— 1|2 



L' 2 (0,t;H) 



+ 



Vp 



H 



E 



j\\L2(0,t;H) 







2 ,t 


+ \\ W \\h(0,t;H) + 


[ Vw 


< c / (^(w 1 


Jo 


H JO 



Then, getting advantage of the Lipschitz continuity of 7 (cf. assumption (13. 24[) ) and 
using (I3.29p . we get 



L 2 (0,t;H) 



+ 



Vp 



+ 



\ W \\l 2 (0,t;H) + 



2 2 

t 2 



2 H + v 



H 



j\\L 2 (0,t;H) 







II 



2 



holding true for all positive constants n and for some c v > 0. Finally, choosing 
T) sufficiently small and applying a standard version of the Gronwall lemma (cf. [TOl 
Lemme A.2]), we get the desired uniqueness of the j3 and w-components of the solution 
to Problem (P) and this concludes the proof of Theorem 13. 101 □ 



5 Proof of Theorem 3.7 



The following section is devoted to the proof of Theorem 13.71 It will be done by 
steps. First we approximate our Problem (P) by a more regular Problem (P) a , 
then (fixed a > 0) we find well-posedness for the approximating problem using a 
fixed-point theorem and then we perform some a-priori estimates (independent of a) 
on its solution. Finally we find existence of solutions of Problem (P) by passing to 
the limit in Problem (P) a as a \ 0. 

In order to perform the following estimates rigorously, we use here the Lipschitz 
continuous Yosida-Moreau approximation dip" = {<p a )' (cf. [TOl Prop. 2.11, p. 39]) of 
dip. 

5.1 The approximating problem 

Let us introduce here the approximating Problem {P) a of our Problem (P) (cf. equa- 
tions (I3.151 - |3~2"0"]) ). First of all, we take a small positive parameter a and we call 7 CT 
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the following Lipschitz continuous approximation of the function 7(10) = exp(w) in 
(13.71) . i.e. the function 

7 »:=( eXpW ifr - 1/a (5.1) 

I (w — 1/cr) exp(l/cr) + exp(l/a) if r > 1/a. 

Moreover let 5 a be the inverse function of 7^ , i.e. 

5 a {Y{r)) =r Vr Gi (5.2) 

and let 7 CT be a primitive of the function 7% i.e. 

7 CT (r) = 1+1 7 CT (s) rfs Vr G R. (5.3) 
Jo 

Then the following properties of 7 CT hold true (cf. also [61 Lemma 5.1]). 
Lemma 5.1. There holds 

7 CT (r) > 7 ff (r), r(5 a )'{j-) > 1 VrGl. 

Now, we can approximate Problem (P) as follows. 

Problem (P) a . Find /3f, ) and (0 CT , £ CT , p 17 ) with the following regularity 
properties: 

w a } $ a G if^O, T; V) n L 2 (0, T; 1/), (5.4) 
ft", ^eH\0,T;H), ufi, eL°°(0,T;V)nL 2 (0,T;W), (5.5) 
r G L 2 (0,T;#), e H\0,T; H) n L°°(0,T;V) H L 2 (0,T;W), (5.6) 



satisfying 



V + Qid t Pl + f^ft 7 + = a.e. in Q T , (5.7) 
+ dtfil + £w a = ft in V' and a.e. in [0, T], (5.8) 

^f!W+r^f?) = f^^ a.,inQ T) (5.9) 



JBfc'J S \£?2 y V 

0" = 7 <7 K), r e ^ CT (/3 CT ), a.e. in Q T , (5.10) 



and such that 



w CT (0)=w , p CT (0)=p a.e. infi, (5.11) 
/3 CT (0) = (/3f (0), /? 2 CT (0)) = /3 a.e. in Q. (5.12) 

Existence and uniqueness of solutions to Problem (P) a directly follow from 
Theorem 13. 101 We proceed now performing a-priori estimates on this solution uniform 
in a in order to pass to the limit in Problem (P) u as a \ 0, recovering a solution 
to Problem (P) . 
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5.2 A priori estimates 



In this subsection we perform a-priori estimates on Problem {P) a uniformly in a, 
which will lead us pass to the limit in Problem (P) a as a \ and recover a solution 
of Problem (P). Hence, let us denote by c all the positive constants (which may also 
differ from line to line) independent of a and depending on the data of the problem. 
Now, let us come to the (uniform in a) estimates on the solution to Problem (P) a . 

First estimate. Test flSUD by p a , (JE2]) by § a = and Q by d t f3 a . Then, 

sum up the result of the first two tests with the sum of the two components of the 
vectorial equation found by the third test. This, after an integration in time on (0, t) , 
with t e [0, T] , thanks to a cancellation of two integrals, leads to 

\\\p°ml + W\\l^ t , H) + [ rKW)+ f I v<r(r)vr 

z Jn Jo Jn 

2 2 

+ Eh^IIW) + ^Eii v ^WIIh+ / ^{P J {t))<\\\Po\\ 2 H + [ 7>o) 

j= i A j= i Jn 1 Jn 



+ f (n^) + u -Y,\\vp-m\ 2 H + [ <p a w+ f [ iwr- 

Jo z =1 Jn Jo Jn 



(5.13) 



Now, following the line of [EJ (5.5)-(5.7), p. 1583], we can deal with the source term 1Z 
recalling (13.111) and using a well-known compactness inequality (cf. pj2, Theorem 16.4]) 
in this way 

I ' R&°< I 11^)11^(0)11^(5)11^(0)^, (5.14) 
o Jn Jo 

[ [ n^^ciiniUoo^ii^) 1 ^!!^^ 

Jo Jan 

< (\\V^) 1/2 \\h(o,t-, H) + Cdm l/2 \\l W) (5-15) 



for all ( > and for some positive depending also on ||IT||£oo(q) (cf. ( 13.101) ) and 
on Q. Now, collecting estimates (I5.14H5~T5"]) . with ( = 1, using Lemma 15.11 with 
assumptions (I3.8H3.101) . we get the inequality: 



\\\p°(t) \\l + 1| w||i W) + 2 1 r(t) + £ jT ds 

2 2 

<c+ [ ||i?( S )||^(n)||^(5)|U 1( n)^ + c||(n 1/2 |li 2( o^)- 
Jo 
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Finally, a standard version of Gronwall lemma gives: 

ibiLw) + nvpii W) + mn 1/2 (t)\\ 2 H + f + J2 \mwiw 

Jo 

.7=1 ^ 



(5.16) 



Second estimate. Testing (15.81) with w a and using (15.161) . it is a standard matter 
to deduce the following bound 

(5.17) 

whence (I5.8P entails 

\\w a \\m(o,T;V) < c, (5.18) 
while, recalling (13.221) . from (I5.16P we also get 

\\r\\w<g T )<c. (5.19) 

Moreover, testing equation (14.61) by |^ cr | _1 ^ 3 ^ CT , we get 

Un\L*r HQT) <c (i = l,2), (5.20) 

and so, by comparison in (I5.9p . we get 

II^A CT |lL5/3( ,T;iy2.5/3(n)) < c (i = 1,2). (5.21) 

Now it remains only to pass to the limit in (l5.7H5.12l) as a \ . This will be the aim 
of the next subsection. 

5.3 Passage to the limit 

As we have just mentioned, we want to conclude the proof of Theorem 13.71 passing to 
the limit in the well-posed (cf. Subsection 15. ip Problem (P) a as a \ using the pre- 
vious uniform (in a) estimates on its solution (cf. Subsection 15.21) and exploiting some 
compactness-monotonicity argument. Let us list before the weak or weak-star conver- 
gence coming directly from the previous estimates and well-known weak-compactness 
results. Note that the following convergences hold only up to a subsequence of a \ 
(let us say ay. \ 0). We denote it again with a only for simplicity of notation. From 
the estimates (I5.16145~2T1) . we deduce that 

w a -> w weakly star in H\0, T; V) n L 2 (0, T; V), (5.22) 

-> d weakly in L 5/3 (Q T ), (5.23) 

%^Pj weakly in H\0,T;H), (j = 1,2), (5.24) 

u/3J -> upj weakly star in L°°(0, T; 1/) fl L 5/3 (0, T; W 2 ' 5/3 (fi)), (j = 1, 2), (5.25) 

p a ^p weakly star in H 1 (0,T;H) nL°°(0,T;V) C]L 2 (0,T;W), (5.26) 

weakly in L 5 / 3 (Q T ), (j = l,2). (5.27) 
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Moreover, employing the Aubin- Lions lemma (cf. [THl P- 58]) and J2TJ Cor. 8, p. 89], 
we also get: 

w a —> w strongly in L 2 (0,T; H), and hence a.e. in Q T , (5.28) 
/3/->/3j strongly in C Q ([0,T];#)nL 2 (0,T;V) if i/ > (j = l,2), (5.29) 
j9 CT ^p strongly in C°([0,T]; H)nL 2 {0,T; V) and hence a.e. in Q T . (5.30) 

Note that (15.231) . ( 15.281) . and relation ( 15 . 101) imply immediately (cf. [18] Lemme 1.3, 
p. 12]) the convergence 

— > $ = 7(w) a.e. in Qt, 
and by Egorov theorem, we also deduce 

iT -»• i? strongly in L q (Q T ) Vg G [1, 5/3). (5.31) 

Finally, it remains to prove the identification of the maximal monotone graph dip, i.e. 

£ G d(p((3) a.e. in Qt, 

with £ = (£i, £2) an d £j (j = 1, 2) that are the weak limits defined in (I5.27p . In order 
to do that we should verify that 

hmsup /V,/3 a )< / T (£,/3). (5.32) 

cr\0 JO JO 

The proof here is splited into two parts. 

Case v > 0. In this case from (15 . 24f45T25|) . using again [ZD, Cor. 8, p. 89], we also 
deduce 

strongly in L r (0,T;L 9 (fi)), Vr G [l,+oo), q G [1,6), (j = 1,2). (5.33) 
^From this convergence and from (15.271) we immediately deduce that (15.321) is verified. 

Case 1/ = 0. In case v — 0, due to assumption ( 13.61) . we also have that f3 is bounded 
in L°°(Qt) ■ In this case we can prove that f3 a is a Cauchy sequence in C°([0, T];H) . 
Indeed one can take the differences of equations ( 15.91) written for two different indices 
a and a' and test it by the difference vector (3 a — (3 a . Take the differences of the 
two equations (15. T|) integrated in time and test it by p a — p a . Summing up the two 
resulting equations and integrating over (0, t), t G (0,T), we get (as a, 0' \ 0): 

\\p a -p^\h { o,t-,H) + \\W - F')m 2 H< [ (r-r'm-p?) 

A JQt 

<2||r-^'|| Ll(Qt) ^0. 

Hence, we have that (3 a — > f3 strongly in C°([0, T]; H) and weakly star in L°°(Qt) 
and so also strongly in L t (Qt) for all r G [1, +00) , which is sufficient in order to prove 
that (15.321) is satisfied. 
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Hence, in both cases, all these convergences with the identifications made above 
make us able to pass to the limit (as a \ or at least for a subsequence of it) in 
Problem (P) CT finding a solution to Problem (P) and concluding in this way the 
proof of Theorem 13.71 Note that the convergences hold only for proper subsequences 
crfc of a tending to because of lack of uniqueness of solutions (cf. also Remark [3.81) . 
This concludes to proof of Theorem 13.71 n 
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